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P ABSTRACT 

^ We present evidence that the electron fraction of the progenitor white dwarf 

1^^ systematically influences the nucleosynthesis of silicon group ejecta in supernova 

^-H type la explosions. In particular, we demonstrate that the abundances generated 

,—1 when quasi nuclear statistical equilibrium conditions apply are preserved during 

W the subsequent freezeout. Analytical relations and post-processing of ID models 
show ^^Si is insensitive to the electron fraction, ^^S has a nearly linear trend with 



Q_i the electron fraction, and ^°Ca has a nearly quadratic trend with electron frac- 

Q tion. We verify that these trends in the nucleosynthesis are reflected in synthetic 

^ spectra. We then explore a new method, based on the quasi nuclear statisti- 

ci cal equilibrium equations, which suggests that measurement of F28Si, Y32s/^28Si; 
^40Ca/^32S! and y54Fc/^28Si in an observed SNe la spectra can be used to con- 

►^ strain the electron fraction in silicon-rich regions. Measuring these four quantities 

CN from the the synthetic spectra generated by our explosion models is sufficient to 

Q^ recover our input Yg to within 6%. Thus, we offer a method to determine the 

"^ electron fraction of the progenitor from observations. 
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Introduction 



Supernovae type la (henceforth SNe la) are thought to be the evolutionary terminus for 
a class of binary stellar systems ([Whelan &: Iben 1973 van den Heuvel et al. 1992 Kahabka 



& van den Heuvel 1997; Parthasarathy et al. 2007 Meng & Yang 2010), the thermonuclear 



incineration of one or more carbon-oxygen white dwarfs (Branch et al. 1995 Wang & Han 



2012), a primary source of iron in galaxies (Matteucci & Greggio 1986; Tang & Wang 2010 



Bulbul et al.||2012" ), accelerators of cosmic rays and sources of kinetic energy in galaxy evolu- 
tion ( Wang|[20lT Powell et al.||2011 ) and useful tools for measuring cosmological parameters 



(PhiUips 1993 Riess et al. 1998 Perlmutter et al. 1999; Kowalski et al. 2008 Wood-Vasey 



eraL][2008| [Hicken et aL][2009| [Riess et al.|[20TT| [Conley et aL][20TTj [Foley fc Kasenl[20TT 



Foley et al.[[2011[ [Sulhvan et al][2011t [Suzuki et aL][2012[ [Silverman et al.[[2012a[|bp . 

The peak luminosity of SNe la is set by the radioactive decay chain ^^Ni — >■ ^^Co — )■ ^^Fe 
(Arnett|1982 Arnett et al.|1985), and the observed photometric correlation between the peak 



luminosity and the timescale over which the light curve decays from its maximum (Phillips 



1993) is understood physically as having both the luminosity and opacity being set by the 



mass of ^^Ni synthesized in the explosion (Arnett 1982 Pinto & Eastman 2000 Mazzali & 



Podsiadlowski 2006 Kasen & Woosley 2007). When corrected for the correlation between 



peak luminosity and light curve decay timescale, the residual dispersion in SNe la optical 
light curves is ~15% (Filippenko 2005), not all of which can be attributed to statistical 



error. This correction removes the dispersion that is attributed primarily from the difference 
^^Ni mass and other physical effects. For example, the residual dispersion in the Hubble 



m 



diagram is reduced by excluding those SNe la with high-velocity ejecta (Foley &: Kasen[[2011 



Foley et al. 2011 Foley 2012). Furthermore, there appears to be a need to correct for the 



host galaxy (Howell 2011) as some properties of the host stellar population are apparently 



imprinted on the explosion. Accounting for such systematic effects potentially allows for 
more accurate determinations of the distance modulus from the observed light curve and 
spectra. This implies identifying and extracting observable physical effects that may create 
a dispersion between the SNe la light curves. 

Over the last decade a number of observational and theoretical studies have sought 
to uncover such systematic effects from variables other than ^^Ni, with the aim of making 
distance measurements more precise and improving our quantitative understanding of the 
progenitor systems. Examples include the empirical correlations between the spectra and 



light curve (Barbon et al. 1990 Branch et aL 1993 Nugent et al.|1995 Blondin & Mezzacappa 
2006J|Bongard et al.|2006| [Branch et al.|2009t|Nordin et al.|[20TT||Blondin et al.|2011t [Foley fc 
Kasen|2011 Foley|2012 ), the dependence of the progenitor metallicity on the peak brightness 
dHoflich et aL]|1998| [Umeda et al.|1999t [Timmes et al.|2003| [Travaglio et al.|2005t [Ropke et al." 
20061 [Ellis et al.| [20081 jGallagher eFall [20081 |Piro fc Bildsten[ [20081 [Chamulak e"raLl[2008 



Badenes et aLl[2008l [Howell et al.1[2009l [Neill et al. [[20091 [Townsley et aLl[2009l [SuUivan et al 



2010 Bravo & Badenes 2011 Foley & Kirshner 2013), asymmetries in the explosion (Howell 



2001 Kasen et al. 2003 2004 Wang et al. 2008 Kasen et al. 2009 Chamulak et al. 2012), 



central density and carbon-oxygen ratio (Hoflich et al. 1998 Dominguez et al. 2001; Ropke 



et al. 2006 Jackson et al. 2010 Krueger et al. 2012), age of the progenitor (Scannapieco & 



Bildsten[[2005| [Mannucci et"aLl[2006l [Sullivan et al.|2006[ |2010| [Krueger et al.|2012D , neutron- 
rich isotope to ^^Ni ratios (Mazzali & Podsiadlowski 2006[), and the opacity of the overlying 



material (Mazzali et al. 2001 Kasen & Woosley 2007). However, we still lack a consensus 



on the progenitor systems as well as how differences in initial conditions create variances in 
the observed properties of SNe la. 

The composition of the white dwarf should have an effect on the nucleosynthesis that 
takes place during the explosion and thus on the isotopic abundances of the final composi- 
tion. For example, most of a main-sequence star's initial metallicity comes from the CNO 
and ^^Fe nuclei inherited from its ambient interstellar medium. The slowest step in the hy- 
drogen burning CNO cycle is proton capture onto ^^N. Consequently, all catalyst nuclei are 
converted to ^^N when hydrogen core burning on the main sequence is completed. During 
helium core burning the reaction sequence ^'^N(a,7)^^F(/3+,z/e)^^0(a,7)^^Ne converts most 
of the ^'^N into ^^Ne. From this point forward, stars have a net positive neutron excess rj, 
defined as ?7 = 1-2{Z) / {A)=l—2Ye where (Z) is the mean atomic number, (A) is the mean 
number of nucleons, and Yg is the electron fraction. Additional burning stages will be driven 
towards producing more neutron-rich elements at the expense of other elements. 

For example, this net positive neutron excess increases the production of neutron-rich 
isotopes such as ^^Fe and ^®Ni instead of radioactive ^^Ni in the regions of the white dwarf 
that produce most of the iron group isotopes during the explosion. This leads to a linear 
correlation between the birth metallicity and the peak brightness (Timmes et al.|2003 Ropke 



et al. 2006; Foley &; Kirshner 2013). However, the range of birth metallicities Zbirth is not 



large enough to account for the full diversity of SNe la peak luminosity ( Gallagher et al.|2005 



2008 Howell et al.[[2009 ) . This has encouraged exploration of other factors that may impact 
neutronization prior to the explosion. Some explorations have focused on the simmering 
phase, where the carbon+oxygen-|-^^Ne white dwarf undergoes an ~1000 yr convective period 
prior to explosion. This convective region is driven by the ^^C+^^C reaction and extends 
outward in mass from from the core to ~1.2 M0. During the simmering phase the sequence 



I2r 



C(p,7) N(e ,z/e) C is primarily responsible for increasing the neutron excess by an amount 



that depends on the total mass of carbon burned prior to the explosion (Piro & Bildsten 



2008 Chamulak et al. 2008). Such studies have established the existence of a "floor" level 



of neutronization that is larger than the neutronization from the birth metallicity when 
Zbirth ~ 2/3 Z0, where Zq is the solar metallicity. Thus, simmering may mask correlations 
between SNe la properties and the birth metallicity. However, for the purpose of this paper, 
what matters is the white dwarf has a well defined Y^ when it explodes, not how the white 
dwarf achieved that Ye distribution. 

On the other hand, observational SNe la surveys exploring the impact of potential 
metallicity effects invariably use the metallicity of the host galaxy as a proxy for the metal- 



licity of the progenitor white dwarf (Ellis et al. 2008 Gallagher et al. 2008 Howell et al. 



2009 Neill et al. 2009 Sullivan et al. 2010). However, it is well established that there is a 
relatively large scatter in stellar iron to hydrogen ratios, A[Fe/H] ~ 0.5 dex, at any given 



age for stars in the Milky Way (Twarog 1980 Edvardsson et al. 1993 Chen et al. 2000 



Feltzing et al. 2001; RoUeston et al. 2000 PediceUi et al. 2009). For example, Feltzing et al. 



(2001) constructed an age-metallicity diagram for 5828 dwarf and sub-dwarf stars from the 



Hipparcos Catalog using evolutionary tracks to derive ages and Stromgren photometry to 
derive metallicities. They conclude that the age-metallicity diagram is well-populated at all 
ages, that old but metal-rich stars exist, and that the scatter in metallicity at any given 
age is larger than the observational uncertainties. Alternatively, by following the chemical 
evolution of homogeneous galaxy models with the evolution of the supernova rates in order 



to evaluate the metallicity distribution function. Bravo & Badenes (2011) find the mean 



metallicity of SNe la and the metallicity of the host galaxy are tightly correlated when both 
metallicities are measured as the CNO abundance. 

If the composition of the white dwarf has an observable effect on the ^^Ni production 
and thus the SNe la light curve, it could have an effect on other elements as well. In this 
paper we present a new direct method to measure the electron fraction Ye in the silicon-rich 
regions for individual SNe la by using the Si, S, Ca, and Fe features from single epoch and 
multi-epoch spectra. Our method follows from the twin facts that ^^Si, '^^S, ^°Ca and ^'^Fe are 
produced in a quasi nuclear statistical equilibrium (henceforth QNSE) environment, and that 
the abundance levels achieved during QNSE do not change during the subsequent freezeout 
as the SNe la expands. Thus, the QNSE abundance levels of these elements are recorded in 
the spectra. Working the situation in reverse, from the observed abundances we can apply 
the QNSE relations to determine the QNSE abundances and thus determine a reasonably 
accurate measure of Ye in the silicon group producing regions. This method is independent of 
any SNe la explosion model and assumes only that the isotopes are synthesized in a QNSE 
state. The motivation behind constraining the electron fraction is to reduce the residual 



dispersion in the Hubble diagram by correcting for a potentially measurable systematic 
effect. Our secondary aim is to provide rigorous nucleosynthesis constraints that can guide 
the modeling of SNe la synthetic spectra. 

In ^ we present our essential results for the trends in the nucleosynthesis, light curves, 
and spectra. In ^we demonstrate the signatures in the nucleosynthesis, hence the spectra, 
are only due to material achieving QNSE conditions. We also present a method for con- 
structing the electron fraction from observations of ^^Si, ^^S, ^°Ca and ^"^Fe in SNe la optical 
spectra. Finally, in Section |4] we discuss the implications of our results. 



2. Essential results for nucleosynthesis, light curves and spectra 

In this section we summarize our main findings on the nucleosynthesis, light curves 
and synthetic spectra near maximum light. We develop a framework which constructs the 
electron fraction in the silicon-rich material that is independent of any specific SNe la explo- 
sion model. We find the abundance of many silicon group (henceforth SiG) and iron group 
(henceforth FeG) elements derive from QNSE conditions. 



We begin with the W7 model (Nomoto et al. 1984 Thielemann et al. 1986 Iwamoto 



et al. 1999) because such models have been extensively analyzed (e.g. Nugent et al. 1997 



Hachinger et al. 2009 Jack et al. 2011). However, as we demonstrate in ^ the abundance 
trends are model independent and depend only on material reaching QNSE conditions. W7 
is a ID explosion model with a parameterized flame speed that captures the stratified ejecta 
observed in "Branch normal" SNe la (Branch et al. 1993). For our purposes, a model is 



a series of temperature and density snapshots from ignition at time=0 s to homologous 
expansion at time=4.1 s. The W7 model assumes a solar ^^Ne mass fraction uniformly 
distributed throughout the white dwarf. Our 'W7-like" models change the assumed value of 
the uniformly distributed ^^Ne mass fraction. Specifically, we control the value of the electron 
fraction Ye by setting the mass fraction of ^^Ne: 1^=0.5 - X{^^lSie)/22=0.5 - Q-X(^^Ne)Q/22, 
where Q is a multiplier on the solar ^^Ne mass fraction X(^^Ne). In our W7-like models. 



22 



Ne represents the metallicity of the white dwarf from both birth and simmering (Timmes 



et al. 2003; Piro & Bildsten 2008; Chamulak et al. 2008; Townsley et al. 2009). 



The nucleosynthesis ejected from these W7-like models is calculated by integrating a 
489 isotope nuclear reaction network ( Timmes||1999 ) over the thermodynamic trajectories of 
each Lagrangian mass shell. This post-processing calculation is not precisely self-consistent 
because the W7 thermodynamic trajectories have the built-in assumption of an energy release 
from carbon+oxygen material complimented with a solar ^^Ne mass fraction. While changes 



to the abundance oof ^^Ne, hence Y^,, shghtly influence the energy generation rate (Hix & 



Thielemann 1996 1999), burned material still reaches QNSE conditions and our analysis 



should still hold. 



2.1. Local and Global Abundance Trends 

First we address whether or not it is justified to use QNSE to recover the abundances 
reflected in the observed spectra of SN la. As Figure [T] shows the local abundances for major 
SiG and FeG elements between mass shells O.TMq and I.SMq at t=1.125 s (solid colored 
lines) and t=4.10 s (dashed black lines). At t=1.25 s the burning front has just passed over 
the 1.28 Mq mass coordinate and most of the ^^Si, ^^S, and ^°Ca synthesized have reached 
their equilibrium abundances. In this mass region, the peak temperatures (3-5x10^ K) and 



peak densities (2-4x10^ g cm ^) ensure QNSE conditions (Nomoto et al. 1984 Thielemann 



et al. 1986). The choice of t=1.125 s is arbitrary and can be replaced by any epoch in any 
model when the material reaches QNSE conditions. At t=4.10 s the explosion has entered 
homologous expansion and synthesis of all the elements has stopped due to the decreasing 
temperature. Complete freezeout does not occur for ^^Fe below O.95M0 at t=1.25 s due 
to residual weak reactions. Figure [T] demonstrates that abundances generated when QNSE 
conditions apply are preserved during the subsequent freezeout. The abundance levels at 
this epoch is reflected in the observed spectra over subsequent days. Therefore, one can 
reverse the process by applying the QNSE equations to recover Ye and the abundances of 
most SiG and FeG elements from the observed spectra. 



Figure |2] shows the flnal mass fractions of the SiG elements ^^Si, ^^S, and ^°Ca between 
O.8M0 and 1.2M0. Two cases for each element are shown, Q=1.0 and 2.0, representing 1.0 
and 2.0 times the solar ^^Ne abundance respectively. Solid lines represent the results from 
post-processed W7-like models and the symbols represent the results from analytic QNSE 
solutions. Figure [s] shows the final mass fractions of FeG elements ^^Fe and ^^Ni over the 
same mass range and white dwarf metallicities as Figure [2J These figures show '^''Ca and ^^Fe 
have the largest systematic changes (up to a factor of two depepending of the mass shell) 
within the silicon-rich region as the electron fraction varies. This suggests that these isotopes 
are more advantageous to examine in SNe la spectra to distinguish between different SNe la 
having different electron fractions. Usually it is difficult to extract the ^^Fe abundance from 
the iron lines. However, ^^Fe is the only iron isotope which is abundant in the regime where 
both ^^Si and ^^S are also abundant in the absence of significant mixing of the QNSE material 
with core material. Therefore, if iron is detected in the ealy time spectra (^ 8 days) at the 
same expansion velocities where SiG elements dominate the spectral features, that iron is 



^^Fe. This result is due to material being in QNSE and is not dependent on any particular 
SNe la model. 

Next we explore the global changes in the abundances of the SiG and FeG elements 
with Ye when our W7-like models reach homologous expansion at time=4.1 s. Fig. |4] shows 
the total ^^Si, ^^S, and '^°Ca molar abundances ejected as a function of the Yg. The curves 
correspond to post-processing the W7-like thermodynamic trajectories and the symbols are 
the results from our analytical QNSE model. Both the post-processing and the model inde- 
pendent QNSE results suggest a nearly constant ^^Si yield with respect to Ye, a systematic 
quasi-linear ^^S yield with respect to Ye, and a more complex trend for the global abundance 
of ^"Ca with Ye. Among the SiG elements, ^°Ca has the largest sensitivity to the electron 
fraction, in agreement with the trends seen in the local abundances. 



2.2. Synthetic spectra and light curves 



We use the PHOENIX radiation transfer code (Hauschildt et al. 1997; Jack et al. 2009 



Jack 2009 Jack et al. 2011) to produce synthetic light curves and spectra from our W7- 



like models. The thermodynamic profiles of the W7-like models end when the explosion 
reaches homologous expansion, about 4 s after ignition. The density, velocity, and abundance 
profiles are then homologously and adiabatically expanded to 5 days after the explosion using 
analytical expressions that account for the local decay of ^^Ni and ^^Co. This is a reasonable 



assumption for SNe la after the initial break out (Arnett 1982). From day 5 onwards, we 



address LTE radiative transfer through the expanding remnant in 0.5 day increments to 
about 21 days to calculate spectra and the light curve. At each of these 0.5 day increments, 
we solve an energy equation that includes the contribution of the adiabatic expansion, the 
energy deposition by 7-rays, and absorption and emission of radiation. As a result we always 
obtain a model atmosphere which is in radiative equilibrium. From the model atmosphere 
we calculate a spectrum and a point in the light curve. 

Figure [5] shows the synthetic light curves for the W7-like models with Q=0.0, 1.0, 2.0, 
and 4.0 representing 0.0. 1.0, 2.0, and 4.0 times the solar ^^Ne abundance, respectively. The 
change in luminosity between these light curves is mostly from the expected linear change in 



^^Ni produced (e.g., Timmes et al. 2003), which decreases with decreasing electron fraction 
Ye. A change of one unit in Q produces a 7% change in the ^^Ni yield and a similar change 
in the peak luminosity. 

Spectral features are a better diagnostic of the electron fraction Yg because abundances 
from different SiG and FeG elements are correlated by the QNSE relations. Synthetic spectra 



at day 15, near peak luminosity, is shown in Figure [6] for the same Fe values used in Figure 
|5j The SiG material has post-explosion homologous expansion velocities ranging between 
9000-13000 km s~^. The continuum optical depth of the SiG material is between 0.2-0.8, 
which was evaluated from the continuum opacity at 5500 A at peak luminosity. Figure |6J 
suggests that the Ca II line at 4250 A changes the most with Ye when compared to S II and 
Si II lines at 5310 A and 5800 A respectively. The Si II line changes the least with variations 
in Ye. These trends, Ca being the most sensitive, S having a near linear dependence, and Si 
the least sensitive, are a reflection of the nucleosynthesis trends shown in Figure |4j 

Changes in the synthetic spectra shown in Figure |6] are due to changes in the abundance 
of a given element and changes in the thermodynamic conditions (such as temperature), 
which themselves are due to the altered elemental abundance. An equivalent width obtained 
from a spectra is a convolution of these abundance and temperature changes. Radiative 
transfer spectral modeling is thus required to determine the absolute abundance of a given 
element from the spectra. However, the ratio of two equivalent widths removes any strong 
temperature dependence. That is, in a region where where all matter has same temperature, 
the ratio of abundances between two elements will be equal to the ratio of their abundances 
times a function that depends weakly on temperature (see Q. Thus, radiative transfer 
modeling of the spectra is not a prerequisite to recovering abundances ratios. Concentrating 
on abundance ratios also makes our approach independent of any underlying assumptions 
about density-temperature-velocity structure of the ejecta. 



3. Constructing Y^ from the abundances 

In this section we establish the QNSE equations which connect the molar abundances 
to Ye. We then discuss how to reconstruct the Ye in the silicon-rich regions from observed 
SNe la spectra. We conclude this section by showing the principle functional dependencies 
of the SiG element abundances on Yl. 



3.1. Basic framework 

We first select a system consisting of the major SiG and FeG elements to trace out 
the most useful equations connecting the individual abundances and their relationship to 
Ye. We choose ^^Si, '^^S, and ^°Ca from the SiG isotopes, and ^^Ni and ^^Fe from the FeG 
isotopes. As figure [T] demonstrate, these are the dominant isotopes under QNSE conditions. 



as 



Conservation of mass and charge can be expressed 

Fn + n + 28^2831 + 321^328 + 40n0Ca + 54F54Fe + SS^SNi 
Fp + 14F28Si + 16^328 + 20y40Ca + 2QY^iYe + 28y58Ni 



1 



(1) 
(2) 



Prom minimization of Helmholtz free energy there follows the fundamental QNSE rela- 



tion (Bodansky et al. 1968 Hix & Thielemann 1996; Meyer et al. 1998 Ihadis 2007) 



Y> 



A,Z 



Y^, 



fip,T)Yt^'-^'-''K'-'' 
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G 



A,Z 



Gk'.tj V Q 



A-A' 



(3) 



where T is the temperature, p is the baryonic mass density, Ga,z is the temperature depen- 
dent partition function, B is the nuclear binding energy, A^^ is the Avogadro constant, ks 
is the Boltzmann constant, and rriu is the atomic mass unit. 

Specification of T, p, Ye, the aggregate molar abundance of the SiG isotopes Vsig, and 
the aggregate molar abundance of the FeG isotopes Ipec is sufficient to solve for all the 
abundances in a two-cluster QNSE environment. There are four independent quantities, as 
IsiG ^-iid Y-pcG are related by equation [3] In this paper, we show that measurement of four 
quantities 1^2881, ^328/^28Si, ^40Ca/>^32S, and Y54Fe/>^28Si from the spectra of a SNe la, is an 
equally sufficient basis to solve for all the abundances in the a silicon-rich region. 



For our choice of isotopes, equation [3] leads to 

^288i 



^328 

^58Ni 
^328 
^OCa 



y-2y-2 
" p n 



lY-'Y-' 



6Y'*Y~ 

" p n 



(4) 
(5) 
(6) 



where a, 7, and S contain the temperature and density dependent terms. 

Next we summarize some simple but essential conclusions which follow from the al- 
gebraic manipulations on the abundances of a few of the major elements. All the molar 
abundances are the local abundances that correspond to a region of the star associated with 
a specific p and T. Using equations |4] and [6| consider the local SiG element ratio 

Y2m(Y^c.Y'^^ ( R 

^328 V ^328 / v5 



<I>(T) 



exp 



knT 



(7) 
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where R = -B28Si ~ -B32S ~ 0.5(i?32s — -B40Ca) and assumes all ratios of nuclear partition 
functions are unity. The quantity $ is independent of the mass density p and contains a 
weak dependence on T. At temperatures typical of QNSE conditions, 2-5x10^ K, fc^T >> R 
leading to $ ~ 1. More precisely, $=0.76 at T=5xl0^ K and $=0.52 at 7=2x10^ K. 
Typical temperatures in the QNSE regions before freeze-out in our W7-like simulations 
are 3.4-4.0x10^ K (a 15% change) inducing a 5.5% change in $. Measuring $ from the 
equivalent width ratios 1^28Si/^32S and l4oca/^32S allows a determination from single epoch 
spectra if the silicon rich material was produced in a QNSE state. With sufficient precision, 
measurement of $ allows a determination of the temperature in the QNSE region before 
freeze-out. We assume for the remainder of this paper that such precision is available and 
the QNSE temperature is a known quantity. Measuring $ from multi-epoch spectra of SNe 
la when silicon features dominate the spectrum then allows trends in the QNSE temperature 
to be assessed. 

A weakly varying $, when treated as a quasi-constant, implies d^/dT ^ 0. It follows 
that the constituent ratios l28Si/^32S and y32s/^40Ca extremize (minimize or maximize or are 
constant throughout) at the same temperature. This property is evident in Fig. |2]- when 
there is a deflection in F28Si/^32S with mass coordinate there is a corresponding variation 
in y32s/5^oca with mass coordinate. In addition, the individual abundances 1^28Si5 ^323, and 
^40Ca extremize at the same (p,T) point where the ratios extremize. In general, isotopes 
with A = 2Z will have a QNSE abundance Iaz that scales as (YpY^)'"', where n is a positive 
or negative integer. These isotopes can be expressed as a function of ^2831/^323 (equation [7] 
is one example) and will extremize at the same (p,T) point as l283i- The only assumption 
in deriving these properties is that the system achieves QNSE conditions, where most of 
the SiG elements are synthesized. The values of the ratios ^28Si/^323 and Kj23/^40Ca at the 
extremum will, of course, depend on the details of the explosion. 



For the FeG elements, equations |4] and |5] imply 



F54Fe ^^323 7 f Rl 

W = — — — = — ~ exp 



y58Ni >^283i a V ksT 

where Ri = (-B323 — -B28Si) — (-B58Ni — -B54Fe)- The quantity \I' is also independent of p 
and contains assumes all ratios of nuclear partition functions are unity. At a typical QNSE 
temperatures, ^(3.4x10^ K)=6.8 and \Ef(4xlO^ K)=5.12 reflecting a 24% change in \Ef for 
a 15% in temperature. Treating \l/ as a quasi-constant and following arguments similar for 
the SiG elements, one flnds y54Fe, y58Ni and F283i also extremize at the same (p,T) point, as 
evidenced in Fig. [3] 
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3.2. A recipe to construct Y^ from the major elements 

For our simplified system consisting of a few major SiG and FeG elements, equation |2] 
may be written as 



Y. = Y, 



28Si 



14 + 16^^ + 20^ + 26^ + 28^ + 



K 



Y, 



28Si 



Y, 



28Si 



Y, 



28Si 



Y 



28Si 



Y, 



28Si 



(9) 



We choose y28Si as the basis because, as the preceding analysis shows, the silicon yield is the 
least sensitive to changes in p, T, and Y^. in the QNSE material. 

The first step in reconstructing Fg is to determine the 1^28Si/^32S and Kioca/^32S abun- 
dance ratios from their equivalent width ratios from a SNe la spectra at the same or nearly 
similar expansion velocities to ensure that the material is coming from the same regime of 
the explosion. The second step is to test if the ratio 

-1/2 

= $ ~ 1. (10) 



Y, 



28Si 



Y. 



328 



Y. 



328 



Y, 



40Ca 



If the ratio is near unity, this verifies the SiG elements were synthesized in a QNSE en- 
vironment. In addition, as described above, this relation may be inverted to determine 
the temperature of the QNSE environment when the SiG elements were synthesized. Mea- 
surement of the ^288i/^32s and y4oca/^32S ratios determines the second and third terms of 
equation [9] since 



K = Y, 



28Si 



-. . , T^^32S , on^40Ca ^328 , or>'^54Fe , oo'^58Ni 
14 + lo-- h 2U— 7 r- h 2d-- h 28 



:iii 



^288i ^328 ^288i ^288i ^288i . 

where we drop Yp since it much smaller {Yp < 10^^) than the other abundances. 

The third step is to measure the y54Fe/^288i abundance ratio. Since ^^Fe is the only 
major isotope of iron synthesized in the same QNSE regions where the SiG elements are 
synthesized in the absence of significant mixing of QNSE material with core material, iron 
measured in an SNe la spectra at the same expansion velocity as l28Si is y54Fe. The fourth 
step is to determine the y58Ni/^288i abundance ratio. From equation |3] one may derive the 
relation 



Y 



58Ni 



Ki 



Y. 



328 



^54Fe ^288i 

exp(-0.0056MeV^/A;BT) 



(12) 
6.7 and 5.0 respectively for 



where the quasi- const ant Ki 
T=3.4xl0^ and 4.0x10^ K. Thus, the last two terms of equation 11 are completely de- 
termined and may be written as 



Y. = Y 



28Si 



-, A , -in ^328 , nn ^40Ca ^328 , or^ ^54Fe , ^o T^ ^328 ^54Fe 

14 + I077 h 20^77 77 h 207^ h 28Ai 



Y, 



288i 



^328 ^288i 



Y 



288i 



^288i ^288i 



(13) 
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Equation 13 is one of our key results. Measurement of four quantities ^2881, ^32s/^28Si) 



^40Ca/^32S) and y54Fe/^28Si ^re Sufficient to determine Y^ from equation 13 for for our sim- 



plified system consisting of a few major SiG and FeG elements. Using this technique and an 
equivalent width analysis on the synthetic spectra shown in Figure [6} we recover the input 
value Yp in our W7-like models to within 



70. 



3.3. Including non-major elements to refine the Yg estimate 

The abundance of any SiG element can, in general, be recovered using equation |3] and 
Y32s/1^28Si, which can be used to improve the accuracy of the Ye determination. For example, 
the abundance ratio of the Ye=0.5 «-chain isotopes y32s/^36Ar is related to ^2881/^328 by 



Y, 



288i 



Y, 



■'328 ■'36Ar 

where K2(T) = exp {0. 0037 MeV/ksT) and at QNSE temperatures is close to unity. The 
values oi K2 are 0.33 and 0.39 for T=3.4xl0^ K and 4.0x10^ K respectively, indicating a 15% 
range in Ki. Including YseAr in the sum for Yg in equation^ adds the term 36/^2(^328/^2881)^ 



to the right hand side of equation 13 Other Ye=0.5 a-chain non-major elements may be 
added in a similar manner. 

Deviation of Y^ from 0.5 comes from the major element y54Fe and additional non-major 
SiG and FeG elements. For a more accurate determination of Yg, one can use mass conser- 
vation, equation [ij to determine the abundances of these non-major isotopes. For example, 
consider the explicit case when FssNi and Ysosi are to be included. Treating Yp and Y^ as 
trace abundance in the QNSE regions, equation [l] becomes 

28r288i + 32F328 + 40y40Ca + 54F54Fe + 5SY,8m + 56F56Ni + SOFsOSi = 1 



or 



Y 



288i 



28 + 32^ + 40^ + 54^ + 58^^^^^ 



+ Y, 



56Ni 



56 + 30 



Y 



30Si 



Y, 



56Ni 



(15) 



(16) 



F288i ^288i ^288i ^288i . 

which we are going to solve for YseNi abundance in the QNSE region. From equation |3] the 
abundance ratio F308i/^56Ni niay be written as 



Y 



30Si 



Y 



K- 



56Ni 



^288i Y^sm _ ^2881 ^58Ni 
Y^mi ^SeNi ^SeNi ^2881 



(17) 



where log(i^3) = -30.61og($). Substituting equation 17 into equation 16, using previously 
developed relations, and rearranging gives 

^328 



56F56Ni + 30 



^28Si ^58Ni 
^56Ni ^28Si 



1-Y 



288i 



28 + 32 



Y 



288i 



13 



+ 40^77 — h 54:^7 h 58Ai 



^32S ^28Si 



Y, 



28Si 



^28Si ^28Si 



;i8) 



Let c = 30 1^28Si ^58Ni/^28Si and b be the right hand side of equation 18 Multiplying equation 
18] by YsGNi thus yields 

fe^GNi + c = (19) 



56^56Ni 



a simple quadratic equation for the VseNi abundance in terms of known quantities. Taking 



the positive root of equation [19] as the correct solution for FseNi, the abundance of Fsosi niay 
then be recovered from equation [17] The derived YsgNi and Fsosi abundances may then be 



used to refine the estimate for Ye by adding the now known terms to equation ^3_ 

Ys2S ^54Fe 



Y, 



Y 



28Si 



14 + 16^ + 20^^ + 26^ + 28K, 



+ Y 



56Ni 



28 + 14 



^28Si 
^30Si 



^328 ^28Si 



Y 



28Si 



^28Si ^28Si 



Y, 



56Ni 



(20) 



Abundances of other non-major elements may be added in a similar manner to the example 
given above. Using this method for the non-major isotopes and an equivalent width analysis 
on the synthetic spectra shown in Figure [6] we recover the input value Yg in our W7-like 
models to within 0.9%. 



3.4. On the functional dependence of abundances on Ye 

In this section we show the trends in the nucleosynthesis that are reflected in the syn- 
thetic spectra — the insensitivity of F28Si to Ye, the near linear dependence of 1^2881 on Ye, 
and near quadratic dependence of Kioca on Ye — follow from the QNSE equations. We begin 
by recasting equation [3] as 



1 



Yaz 



yA-a'-{z~z')yZ~z' 
f{p,T)YA'z' " " ■ 

We may assume without loss of generality that A > A' and Z > Z'. Now let 

iT^protons 



(21) 



(22) 



which is identical to the definition of the electron fraction Ye but without the free protons. 

(23) 



With 6=YA'z'fip,T), equation 21 may then be written 

Ya^Z, = 9{1 -Ye- wf^-^'iYe - Wf^-^' 
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Multiplying this equation by Zj and summing, 



w 



i^protons i^protons 



(24) 



The right hand side of equation 24 is significant only when Zi is closest to Z' . For the SiG 
group pair ^^Si and ^^S, equation 24 becomes 



w = eZ^2s{l-Ye-wf{Y,-wf 



which is a quartic in w. To order Y^, equation 25 has the solution 



w 



2ez- 



325 



(2y, + 1) ± ^/{AY,ez,2s + 1) 



(25) 



(26) 



Substituting Z'=16 and expanding the square root in equation [26] leads to a zeroth order 
term that is a constant independent of Y^ and a first order term that is is linear in Y^. From 



the expression for w in equation |24[ the largest contribution comes from 2'28Si5^28Si- We thus 
identify F28Si as the constant term and the linear term in Y^ with 1^323- Finally we identify 



^40Ca with keeping the higher order terms in Y^ in equation 25 



Discussion 



We have presented several lines of evidence that the electron fraction systematically 
influences the nucleosynthesis of the SiG elements in SNe la. In particular, we demonstrated 
that the abundances generated when QNSE conditions apply are preserved during the subse- 
quent freezeout. We showed the explosion model independent QNSE relations predict ^^Si is 
insensitive to Y^, ^^S has a near linear trend with Yf,, and the ^*^Ca yield has a near quadratic 
trend with Yg. We verified these trends in the nucleosynthesis are reflected in the synthetic 
spectra of W7-hke models. Finally, we suggested that measurement of 1^2881 and the abun- 
dance rations ^323/^2881; ^40Ca/^32S) and V54Fe/^28Si in an observed SNe la spectra can be 
used to constrain the electron fraction in the silicon-rich region. Our technique allowed us 
to recover the input Y^ of our W7-like models to 6%. Including other non-major SiG and 
FeG elements allowed us to recover the input Ye within 0.9%. The reconstruction of Ye in 
the silicon dominated layers is not based upon assumptions of adiabatic conditions or the 
idiosyncrasies of a particular numerical simulation. We emphasize that we use our W7-like 
models to verify the claims but the claims are not derived from such W7-like models. 

Finally to address concerns about how the abundances be obtained from the equivalent 
widths including ionization, we outline the relevant relation. Let Ntot be the total number 
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of atoms for a given element or its abundance, Nj be the number of neutral atoms of the 
given element, Njj be the singly- ionized atoms of the same element, N^ be the number of 
singly-ionized atoms in ground state responsible for the equivalent width corresponding the 
given transition, and A^^ be the number of singly-ionized atoms of the same element which 
are in the excited state and hence do not contribute to the absorption profile. The ratio of 
excited to ground state population (for the singly-ionized) atoms, br, is calculated from the 
Boltzmann equation. The ratio of singly-ionized to neutral atoms of the same element, r^, is 
calculated from the Saha equation. If an absorption profile is saturated then the equivalent 
width EW (X foscNd where fosc indicates the oscillator strength of a given line transition. 
For a resonant transition fosc = 1- We have 

Ntot = Ni + Nii 



NhH + 



-1^ 



= (iV„ + iV,)(l + r, 

= N,{br + l){r:' + l) (27) 

The equivalent width of Si, for example, is oc foscNd,28Si where fosc is the oscillator strength 
of the observed line transition and iVd,28Si is the number of Si atoms contributing to Si II 
line. It follows that 



{br + l)(r;i + l)EWsi (X Xsi * /si = 28Fsiexp (-^] , 



(2^ 



where /si is the ionization fraction for singly ionized silicon, and e denotes the ionization 
energy. To be more precise, this expression should also have multiplicative Gaunt factors 
involving the cross section of a given line, which for resonant lines are very close to unity. 

The key relation for determining if QNSE conditions apply, equation [7} in terms of 
equivalent widths becomes 

EWs^fEWc^Y'' _Xs^fXcA'^'' f e{St)-e{S)-Q.h[e{Ca)-e{S) \ 

^-' EW, \EWJ ~ X,\Xs) '""P I ksT 

(29) 

where /corr = /s,corr /r,corr- The quantity /r,corr gives the correction due to non-zero level 

population of higher ionization states and fs,corr gives the correction due to ionization. In 

terms of known quantities, and 



1 + 6r,28Si / 1 + ^r,40Ca 
Jr, corr 



1 + &r,32S V 1 + ^r,32S 

1 

1 + lAs,28Si /I + l/^s,40Ca^ ' 



^''''°" " 1 + l/r.,32S V 1 + 1/^.,32S ' ^ ^ 
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Inserting numerical values for the ionization energies, the exponential term in equation 27 
is exp{0. 087 /{KbT)). For T=7000 K, the exponential is ^ 1.15, while for earlier epochs at 
T=10000 K, the exponential is ~ 1.1. Non-LTE effects and line blending can potentially 
impact these relations between the equivalent width and element abundance level. 
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Fig. 1. — Mass fractions of the major elements during a W7-like explosion. Solid colored lines 
correspond to t=1.25 s and and dashed black lines correspond to t=4.10 s since ignition. The 
abundances produced by QNSE conditions at t=1.25 s are the same as the final abundances 
when nuclear reactions freezeout at t=4.0 s. The expansion velocities on the upper x-axi are 
when the W7-like model reaches peak luminosity. 
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Fig. 2. — Mass fractions of the major SiG elements as a function of interior mass for the 
post-processed W7-hke models (solid and dashed lines) and the QNSE calculations (symbols) 
at t=1.19 s. Two cases are shown, one for a solar ^^Ne abundance {Q=l) and one a twice 
solar ^^Ne abundance {Q=2). Overall, the agreement between the post-processed W7-like 
models and the analytical QNSE results are satisfactory. '^°Ca shows the largest change, up 
to a factor of two, while ^^Si is insensitive to changes in ^^Ne. 
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Fig. 3. — Same as Figure |2] but for the major FeG elements. ^^Fe shows the largest change 
with changes in ^^Ne. Note that ^^Fe is the only significant iron isotope present when the 
SiG elements of Figure |2] are dominant. 
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Fig. 4. — Global abundances of ^^Si, ^^S, and ^°Ca as a function of the electron fraction 
Ye produced by the post-processed W7-like models (hues) and the analytical QNSE results 
(symbols). As for the trends in the local abundances (Figs. [l]-|3]), ^^Si is independent of Yg, 
^^S shows a near linear dependence, and '^°Ca shows a more complex, but near quadratic, 
dependence with Y^. 
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Fig. 5. — Light curves for the W7-hke models with 0.0. 1.0, 2.0, and 4.0 times the solar 
^^Ne abundance. 
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Fig. 6. — Spectra for the W7-like models with 0.0. 1.0, 2.0, and 4.0 times the solar ^^Ne 
abundance. The Sill, SII, Fell, Mgll and Call H+K lines at rest frame wavelengths 5800, 
5310, 4750, 4610, 4250 A are labeled. 
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